Spectral behavior of contractive noise 
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We study the behavior of the spectra corresponding to quantum systems subjected to a contractive 
noise, i.e. the environment reduces the accessible phase space of the system, but the total probability 
is conserved. We find that the number of long lived resonances grows as a power law in h but 
surprisingly there is no relationship between the exponent of this power law and the fractal dimension 
of the corresponding classical attractor. This is in disagreement with the predictions of the fractal 
Weyl law which has been established for open systems where the probability is lost under the effect 
of a projective noise. 
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I. INTRODUCTION 

In recent years there has been a great interest in the 
study of open quantum systems. The motivation for this 
upsurge can be found in the recent development of quan- 
tum computation and information 0, [J (even though 
some related aspects have been under investigation from 
long ago, like in scattering theory, for example). The 
noise coming from the environment has always been the 
major drawback for any realistic attempt to implement 
a quantum computer and a serious obstacle for quan- 
tum information purposes. At the same time we have 
realized that we know much less about quantum open 
systems than about their closed counterparts. There- 
fore, their study has become an extremely active topic in 
fundamental physics Q. 

It is in this context and for the case of scattering 
systems that the fractal Weyl law for the number of 
long-lived resonances has been conjectured and tested by 
means of several examples The classical invariant 
distribution in these cases, i.e. the fractal hyperbolic set 
of all the trajectories non-escaping in the past and fu- 
ture (the repeller), plays a fundamental role with respect 
to the quantum spectrum of quasibound states (the res- 
onances). Indeed, the aforementioned law says that [5[ 
the number of resonances that decay at the slowest pace 
(the long-lived ones) grows as h~( 1+dH \ where dn is the 
partial Hausdorff dimension of the repeller. We can trace 
back these studies to the proposal and proof of a fractal 
Weyl upper bound for a Hamiltonian flow showing a frac- 
tal trapped set @. Later, results that seem to strongly 
support the validity of the conjectured law have been ob- 
tained for different kind of systems ranging from smooth 
to hard walls potentials Q- But it is in open quantum 
maps that the fractal Weyl law was more easily tested [8| . 
In these systems, paradigmatic in the quantum chaos lit- 
erature, the resonances have been found to grow as hr d , 
where d is the partial fractal dimension of the repeller. 
However, explorations of the spectral behavior of quan- 
tum systems having different kinds of invariant classical 
distributions associated with them are very scarce. In 
this sense it is very interesting to ask what happens in 
the case of a non-pro jective kind of noise, i.e., a contrac- 



tive environment, for example. Here we will focus on dis- 
sipative quantum operations, whose action can be inter- 
preted as aphase space contraction leading to dissipative 
dynamics [9] . There were very few attempts to study the 
spectral properties of these kind of systems in previous 
works. The only antecedent we were able to find in the 
literature is the paper by Ermann et al. [Toj . But in this 
case only the spectrum of a discretized Perron-Frobenius 
operator was considered. For this approximation to the 
classical problem (i.e. not a proper quantization of it), 
the authors found that the behavior of the corresponding 
long-lived resonances follow the fractal Weyl law. 

In this work we analyze the dissipative baker map. To 
obtain the quantum counterpart we have implemented a 
standard procedure in which the noise superoperator is 
written in terms of appropriately defined Kraus opera- 
tors [ll|. For this model, all classical initial conditions 
fall on a strange attractor after a few iterations. At the 
quantum level, this is represented by a resonance with 
eigenvalue 1, and the rest of the spectrum lying inside 
the unit circle. In sharp contrast to what has been found 
for the discretized classical dynamics, the number of long- 
lived resonances of the superoperator grows as a power 
law in H, but the exponent is rather insensitive to the 
dimension of the fractal invariant set. 

This paper is organized as follows: in Section |H] We 
present the definition and details of the model we have 
used to study this kind of spectrum. In Section IIIII the 
numerical results are analyzed and we present an inter- 
pretation of them supported by further studies. Finally, 
we draw the conclusions in Section IIVI 



II. MODEL SYSTEM 

One of the most simple models one can think of at the 
time of studying complex systems are chaotic maps. Re- 
gardless of their simplicity they capture all the essential 
features of chaotic behavior and their quantization allows 
to extend these advantages to quantum mechanics. All 
this turned them into paradigmatic models for quantum 
chaos and the theory of dissipative systems |L2h15| . We 
have investigated the spectral behavior of the dissipative 
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baker map. The classical map is defined on the 2-torus 
T 2 = [0, 1) x [0, 1) by 



B(q,p) = 



_ f (2g,ep/2) 



(2 g -l,(ep+l)/2) 



if < q < 1/2 
if 1/2 < g < 1 



(1) 



This transformation is an area-contracting, piecewise- 
linear map. The map contracts the torus in the p direc- 
tion by a e factor, stretches the unit square by a factor 
of two in the q direction, squeezes it by the same factor 
in the p direction, and then stacks the right half onto the 
left one. After a few steps any initial distribution falls 
into a fractal set, the strange attractor. 

In an even iV-dimensional Hilbert space (where N — 
(27rfi.) _1 ), the quantum baker map is defined in terms of 
the discrete Fourier transform in the momentum repre- 
sentation as [H, [13] 
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(2) 



with 



(G N ) jk = -= exp{-27ri(j + l/2)(fc + l/2)/N}. 
V N 



Bn is a unitary matrix and represents the quantum dy- 
namics of the closed baker map. 

As discussed in [9| quantum dissipative processes can 
be described by non-unital quantum operations. In this 
work the dissipative noise is implemented by an N 2 x TV 2 
Kraus superoperator of the form: 



III. RESULTS 

The spectrum obtained by the diagonalization of the 
superoperator of Eq.© consists of a leading eigenvalue 
A = 1 , and TV 2 — 1 complex eigenvalues A (with 
|A| = exp(^)) inside the unit circle. Due to its non- 
normality, $ has distinct left ip^ and right eigenoperators 
ipjf corresponding to each eigenvalue. As expected, the 
Husimi function of the invariant state i^o having eigen- 
value A = 1, closely follows the fractal structure of the 
attractor corresponding to the classical dissipative map. 
This is shown in Fig. [T] for several values of e. In the 
classical case the dissipation parameter e determines the 
dimension of the attractor d = 1 + (ln(2)/(ln(2) — ln(e)) 
[l8j . In the quantum case e is related to the quantum 
phase space contraction rate @ . We have verified that as 
N increases, the phase space representation of i/jq reveals 
finer details of the attractor , reflecting the quantum to 
classical correspondence. 



JV-l 

M = ^ A * ® A ^ (3) 

M=0 



where 



A " = E \ ( i - ^ ) - e )*|i - A*><*| (4) 

are operators accounting for transitions towards the mo- 
mentum state \i = 0). The coupling constant e coincides 
with the dissipation parameter of the corresponding clas- 
sical map. M is a trace preserving QZ M -A^-A^ = 1) 
and non-unital A^A^ ^ 1) superoperator, which de- 
scribes a process in which phase space volume is not pre- 
served. Superoperators of this type are obtained from the 
integration of master equations derived from modelling 
a microscopic interaction of a system (an oscillator, a 
large spin, etc) with a thermal bath at zero temperature 
representing the environment [l9l - |2~i1 |. 

To describe the noisy evolution of the density matrix 
of the system we compose M with the unitary map ©, 
according to: 

$ = (B N <8> Blf) o M. (5) 
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FIG. 1: Husimi representation of the invariant state of the 
quantum dissipative baker map for (from top to bottom) e = 
0.8, 0.6, 0.4 and N = 180 (left column), and the corresponding 
classical attractors (right column). 



As a first approach to analyzing the behavior of the res- 
onances, in Fig. [5]we show different spectra of $ for three 
values of the dissipation parameter e = 0.8,0.6,0.4 and 
a Hilbert space dimension N — 180. We notice that the 
longest-lived eigenvalues, though change their position in 
the complex plane, more or less keep their moduli while 
the radius r\ of the dense circle where most of the eigen- 
values concentrate strongly shrinks for increasing values 
of the dissipation. However, this radius cannot be di- 
rectly related to the parameter e, in contrast to what is 
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observed in [T(| for the spectra of the discretized Perron- 
Frobenius operator. This gives to the quantum version a 
seemingly more contractive character in spectral terms. 
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FIG. 2: Eigenvalues of the quantum dissipative baker map $ in 
the complex plane for three different values of the dissipation. 
From left to right: e = 0.8, 0.6, 0.4. In all cases N = 180. 
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To better show this point the differential radial density 
distributions defined as A ^ 7 ^ are plotted in Fig. O for 
different values of N. These densities, normalized accord- 
ing to Jp 18 A ^ 7 ' 1 = 1, are practically N-independent for 
the dimensions we are considering. They present a rela- 
tive maximum at a decay rate 7 -C — 2 In e, and then they 
smoothly increase at larger values of 7. These profiles 
are very different from the ones obtained in [13] for the 
spectrum of the Perron-Frobenius operator, which peak 
around 7 = — 21ne, that is, the global relaxation rate to 
the strange attractor, and decay rapidly for larger 7. 
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3: Density of states A ^ 7 ^ as a function of 7 for different 
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values of N (90(a), 100(D), 150(«), 180(*))and two values of 
the dissipation e = 0.8 left and e = 0.6 right. The densities 
are normalized according to J Q 18 A ^ 7 ^ = 1. 



We now come to the central point of this paper. We 
want to see how the differences between a full quantiza- 
tion and a discretization procedure of phase space which 
are suggested by the comparison of the corresponding 
spectra affect the adherence to the fractal Weyl law. Fig. 
S] displays the fraction of long-lived resonances Nl< ^'i ut 
as a function of N 2 for e = 0.8 and values of "f cut ranging 
from 4 to 20. It can be clearly seen that this fraction 
scales as jV ' y ^l c " t oc (N 2 )~P with an exponent ft ~ 0.8 
fairly independent of the cut-off value 7 cu t used to dis- 
tinguish between the long-lived states and the short- 
lived ones. The existence of a power law dependence 



FIG. 4: Log-log plot for the fraction of long-lived resonances 
(with 7 < 7cut) as a function of N 2 for e = 0.8 and eight 
values of the cut-off value ranging from 4 < j cut < 20. The 
fitted slopes are in the interval 0.76 < ft < 0.82. 



for N ~' < ^"i ut and the insensitivity of the exponent to the 
choice of the cut-off value 7 c „t are fundamental properties 
of the fractal Weyl law which has been originally conjec- 
tured for projective noises. However, in the case of our 
contractive noise the exponent which fits the numerical 
results cannot be related in a direct way to the dimen- 
sion of the classical attractor. This dimension is equal to 
d = 1.756 for e = 0.8 and thus the relation ft — 2 — d 
does not hold. 

Moreover, the exponent ft turns out to be practically 
insensitive to the value of the dissipation parameter e. 
This remarkable fact can be appreciated in Fig. [5] where 
JV ' T ^2 c " t (N 2 ) is shown for different values of the coupling 
constant e = 0.8,0.6,0.4 (and three choices of jcut)- In 
all the cases the data follow a power law, with ft ~ 0.8, 
independently of the fractal dimension of the correspond- 
ing classical attractor ( d — 1.756 for e = 0.8, d = 1.576 
for e = 0.6, d = 1.431 for e = 0.4). 

For comparison we have applied the superoperator for- 
malism to the well-studied case of an open Baker map 
subjected to a projective noise modelled by a non-unitary 
operator U open , by defining $ ope „ = U op en ® U% pen . We 
verified that, as expected, the fraction of long-lived res- 
onances of $ op en scales with iV 2 as does the fraction of 
resonances of U open with N (we have taken N values in 
the same range than for the contractive case), thus fol- 
lowing the fractal Weyl law. 

In order to shed more light on these results we turn 
to analyze the eigenvectors. To characterize the eigen- 
states of $ we consider the overlap of the right eigen- 
states with the invariant state, by defining the measure 
KV^o 3 )! = \Tr{^^)\. Fig. M displays the depen- 
dence of this measure with respect to 7, for different 
values of e. It is clear that the overlap is in the aver- 
age larger for slow-decaying states, and decreases as the 
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FIG. 5: Log-log plot for the fraction of long-lived res- 
onances (with 7 < 7cut) as a function of iV 2 for e — P 
0.8(*),0.6(»), 0.4(D), with a) j cut = 20, b) j cut = 15.2, c) 
■y cut = 9.2 . The fitted slopes are /3 = 0.83,0.75,0.78 for a), 
/3 = 0.87, 0.76, 0.82 for b) and p = 0.84, 0.82, 0.77 for c). . 
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decay rate increases. But it is always small, this giving 
further support to the strong contraction of the spectrum 
we have observed. The distinction between slow- and 
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0.125'- - 
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FIG. 6: Overlap \(ip^if>o')\ as a function of 7 for different 
values of the dissipation parameter e = 0.8, 0.6, 0.4 (from left 
to right) and N = 180. 



FIG. 7: Husimi representations of the a) left and b) right 
eigenstates for the closest to the unit circle excited state |A| = 
0.5845 (upper panel) and for a state with a large decay rate 
|A| = 9.278 x 10" 5 (lower panel), e = 0.8, N = 180. 











* * m » m « - 




- • mm .5 




v - r 



q 1 



fast-decaying states is also noticeable if we look at their 
distribution in phase space. For this, we compute the 
Husimi representation of the right and left eigenstates of 



$, (z\^ >Ij \z) = Tr(ip"' \ \z)(z\) where \z) are coherent 
states centered at z — (q,p). As shown in Fig. [7]for a typ- 
ical slow-decaying state (upper panel) the Husimi density 
of the right eigenstate reasonably follows the structure 
of the attractor (though being much more localized than 
the invariant state), while the left one has a delocalized 
pattern. For a state with a large decay rate (lower panel) , 
the probability pattern becomes difficult to associate to 
the attractor. 

Although the long-lived eigenfunctions are morpholog- 
ically different from the invariant state, we have verified 
that they have support on the phase space region cor- 
responding to the classical attractor by making use of 
a recently developed representation especially suited for 
open systems [22|. We have built the sum 



E 



A=A C 



(zW*\z)(z\ri\z) 



(6) 



and obtained the distribution shown in Fig. |8l indicating 
a clear localization of the long-lived states on the attrac- 
tor region in the semiclassical limit. 



FIG. 8: Husimi representation for the projector of Eq. ([6]), 
with \ cut = 0.1348 (or j cu t = 4) , e = 0.8, N = 180. 



In view of these results we have found that though the 
classical support of the long-lived resonances is the clas- 
sical attractor, their number scales with TV at a different 
pace than that predicted by the fractal Weyl law. 



IV. CONCLUSIONS 

We have studied the spectra of a paradigmatic model 
of the quantum chaos and dissipative systems theories, 
the baker map with dissipation, by following a standard 
quantization procedure based on the Kraus representa- 
tion of superoperators. In contrast to what happ ens 
for the discretized Perron- Frobenius operator in [lOj, we 
have found that the standard fractal Weyl law does not 
hold in this case. Even if the fraction of long-lived reso- 
nances does scale with the dimension following a power 
law which is roughly insensitive to the cut off value of 
the decay rate considered for the statistics, the exponent 
(which is approximately constant) is not directly related 
to the fractal dimension of the classical attractor. 
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In order to give an interpretation for this intriguing 
result we have investigated the morphology of the eigen- 
states of the system. In particular we have built the 
Husimi representation of the projector constituted by the 
eigenfunctions with slow escape rate [22| . We found that 
its density concentrates on the attractor region, strongly 
suggesting that one should expect, as for the projective 
case, a connection between the statistics of these long- 
lived resonances and the fractal dimension of the strange 
attractor. 

A possible cause for the lack of such natural connection 
might be in the non-orthogonality of the eigenfunctions 
localized on the strange attractor. In fact, the fractal 
Weyl law supposes the quasi-orthogonality of the eigen- 
functions supported by the fractal set, since it is based 
on a Planck-cell partitioning of this set. This works in 



the case of projective noises, but might not be the case 
for contractive ones. This is an open problem and we are 
currently studying [23[ the degree of non-orthogonality 
induced by a contractive dynamics as an explanation for 
the apparent non validity of the standard Weyl law re- 
vealed by the present numerical investigation. 
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